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Previous  experiments  on  plasma-propellant  interaction  have  shown  enhanced  bum  rates.  The  temperature  of 
the  propellant’s  bed  may  play  a  role  in  the  bum  rates  under  plasma  injection.  In  this  report,  a  set  of  experiments 
were  first  conducted  on  the  plasma-flow-field  to  determine  the  spatial  and  temporal  distribution  of  the  plasma 
temperature,  pressure,  number  density,  and  velocity.  The  experiments  revealed  a  decreasing  plasma  pressure, 
plasma  temperature  and  plasma  number  density  as  plasma  is  leaving  the  capillary  source  and  expands  in  air.  The 
plasma  jet  velocity  2  inches  from  the  source  exit  was  found  to  be  about  1300  m/s.  Following  characterization  of 
the  plasma-flow-field,  a  set  of  experiments  were  conducted  on  JA-2  solid  propellant  with  controlled  bed 
temperature.  Experiments,  with  proposed  models,  revealed  a  functional  form  that  includes  the  bed  temperature. 
Increased  bum  rates  were  observed  with  increased  bed  temperature,  which  appears  to  follow  a  power  law.  A 
model  for  the  bum  rate  is  proposed,  which  includes  the  bed  temperature  and  has  the  BR  =  A  Pa  (  T/  T  ambient )  b. 
Further  investigation  on  plasma  parameters  helped  in  modifying  the  model  to  include  plasma  radiative 
heat  flux  and  plasma  number  density  in  the  proposed  model. 

I.  Introduction 


Electrothermal  plasma  has  been  studied  as  a  viable  replacement  for  solid  propellant  ignition.  The 
studies  were  not  only  conducted  for  possible  ignition  using  plasmas,  but  also  as  means  of  controlling 
and  enhancing  bum  rates  of  propellants.  With  plasma  augmenting  the  propellant’s  bum  rate,  a  mix 
between  plasma  electrical  energy  and  chemical  energy  of  the  propellant  may  provide  a  possible  means  of 
achieving  high  enthalpy  flows.  Such  high-enthalpy  flows  may  provide  better  interior  ballistics  and 
higher  muzzle  velocities  for  hyper- velocity  launchers  [1-8].  Previous  work  has  shown  that  the  plasma 
ignited  propellant  bum  increases  surface  erosion  compared  to  that  of  conventional  ignition  [1-5].  Those 
experiments  studied  the  effects  of  plasma  pressure  impact  on  the  bum  rate  of  propellant  when  the 
propellant  was  at  ambient  conditions.  No  propellant  bed  temperature  sensitivity  has  been  examined,  as 
of  yet,  when  the  propellant  is  ignited  by  electrothermal  plasmas. 


This  work  will  determine  if  any  propellant  temperature  dependence  exists  in  the  bum  rate  and,  if  it 
exists,  the  functional  form  of  that  dependence.  A  series  of  experiments  was  performed  on  NCSU’s 
PIPE  [2,4,5,  9,10]  ETC  device  that  impacted  plasma  onto  a  heated  propellant  bed  (JA2).  Experiments 
were  conducted  using  an  incomplete  “extinguished”  bum  technique  to  provide  an  evaluation  of  the 
erosive  bum.  These  plasma  impacts  resulted  in  propellant  mass  loss  that  could  be  correlated  to  bum 
rate.  If  propellant  bed  temperature  significantly  impacts  the  bum  rate  of  the  propellant  then  this  would 
be  shown  in  the  experimental  data  returned  from  the  mass  loss  as  compared  to  bum  rates  at  ambient 
temperature.  The  experiments  also  required  a  pressure  difference  so  the  propellant  was  placed  at  1,  2 
and  3  inches  from  the  barrel.  This  allowed  for  a  correlation  based  on  pressure  and  propellant  bed 
temperature.  Additionally,  the  radiative  heating  effect  of  the  plasma  on  the  propellant  bed  could  be 
evaluated  from  the  plasma  temperature  change  farther  away  from  the  barrel. 

The  main  goal  of  the  work  is  to  determine  a  functional  equation  that  can  be  used  with  common 
diagnostics  to  determine  the  impact  of  the  bed’s  temperature  on  bum  rates.  Pressure  is  commonly 
measured  in  expansion  chambers  and  is  relatively  easy  to  obtain.  As  a  result,  the  functional  form  will 
have  pressure  and  bed  temperature  dependence  and  the  bum  rate  with  these  dependencies  can  be 
calculated.  The  final  goal  of  all  current  propellant  studies  is  the  development  of  a  plasma-sensitive 
propellant  that  will  only  combust  when  burned  with  plasma.  This  will  be  a  safer  propellant  that  will  not 
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combust  under  fire  or  by  any  other  accidental  means.  This  work  continues  the  efforts  to  develop  such  a 
plasma  sensitive  propellant  for  the  viable  use  of  ETC  launchers. 

Many  researchers  have  studied,  in  detail,  the  interaction  of  plasma  with  propellants.  Electrothermal 
plasmas  are  used  because  they  are  characterized  by  their  high-density  and  relatively  low-temperature  (1- 
lOeV)  such  that  ET  plasma  resembles  a  blackbody  [7-16]. 

The  majority  of  the  research  at  the  ARL  has  been  performed  in  order  to  determine  the  effects  of 
plasma  temperature,  pressure,  and  radiation  effect  [16-21].  In  general,  these  are  compared  with  the 
conventional  black  powder  ignition  means  in  a  large  assembly  that  simulates  how  a  full  propellant 
packed  charge  would  bum.  These  provide  a  good  indication  of  how  a  real  cannon  device  would  be 
operated  and  an  idea  of  the  geometry  of  the  bum  propagating  through.  However,  these  closed  bomb 
devices  don’t  lend  themselves  well  to  the  study  of  optimization  because  of  the  burning  and  pressure  due 
to  other  grains.  Individual  grain  study  is  useful  for  determining  the  actual  effect  of  plasma  on  the 
propellant  and  gives  an  idea  as  to  the  improvements  in  packing  geometry  or  even  chemical  composition 
to  give  the  desired  bum  characteristics.  NCSU’s  plasma  engineering  group  has  ETC  devices  that  are 
ideally  suited  for  single  propellant  grain  plasma  experiments  and  as  a  result  are  contracted  by  ARL  to 
perform  propellant  studies. 

Previous  work  has  shown  geometry  effects  by  rotating  the  propellant  to  change  the  plasma  impact 
angle.  The  study  revealed  that  normal  impact  results  in  the  highest  bum  rate. 

The  traditional  form  of  the  combustion  equation  comes  from  chemical  kinetics.  For  example  a 
bimolecular  reaction  of  the  type 

A  +  B-+M  +  N 


The  rate  of  reaction  is  proportional  to  the  probability  of  collision  in  the  molecules.  This  probability  is 
determined  by  the  number  of  moles  per  unit  volume  or  the  concentrations.  Hence 

d(VCA)_  dNB  _  d(yCg) 
dt  dt  dt  dt  A  B 


where  V  is  the  volume  and  k  is  the  reaction  rate  constant. 

Performing  differentiation  and  dividing  through  by  a  constant  volume  yields 


=  -^L  =  jfeC  C, 
dt  dt  A  B 


If  the  reaction  proceeds  as  the  stoichiometric  equation  states 

aA+bB  — »  mM  +  nN 


then  the  rate  of  change  of  the  concentration  is  related  to  those  remaining  and  can  be  written  in  the  form 

[22] 
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-kCn 

-  Kn^A 


The  concentration  can  be  written  as  density  times  the  mole  fraction  and  taking  the  full  derivative  of 
the  concentration  yields 

-k(psy 

ot  ot 

Dividing  by  a  constant  density  and  solving  for  the  change  in  mole  fraction  gives 

The  mole  fraction  lost  due  to  combustion  is  proportional  to  the  bum  depth  of  the  material  and 
density  times  the  mole  fraction  is  proportional  to  pressure.  This  allows  the  equation  to  be  written  in  the 
conventional  form  of  the  Saint  Roberts  Equation  [5] : 


BR  =  bPa 


with  P  in  psi  where  ct  takes  care  of  the  combination  of  n  and  n-1.  The  value  for  the  exponent  n  was 
determined  experimentally  for  the  modified  double  base  propellant  JA2  to  be  n=1.3755.  The  b  was 
determined  to  be  8.865x1  O'5.  This  value  is  higher  for  plasma  impacted  propellant  than  for  the 
conventional  black  powder  combustion  n=0.889  and  b=203.2xl0'5  [4],  This  results  in  higher  bum  rates 
for  an  equivalent  pressure  (because  the  value  of  b  is  constant  for  each  propellant)  which  would  change 
the  ballistics  of  a  projectile  launched  by  the  propellant  (and  is  necessary  to  know  in  the  case  of 
cannons).  Bourham  also  showed  experimentally  that  the  bum  rate  is  not  dependent  solely  on  pressure 
but  is  also  coupled  with  the  radiation  heat  flux  from  the  plasma  temperature.  This  was  done  by 
comparing  effects  of  single  and  double  exposed  samples  [4].  The  samples  that  were  exposed  twice, 
once  on  each  face  of  the  propellant,  showed  a  higher  gasification  rate.  The  conclusion  was  drawn  that 
this  gasification  rate  was  due  to  the  in-depth  reactions  caused  by  the  radiation  flux  of  the  plasma. 

It  has  also  been  shown  that  the  only  augmentation  the  plasma  gives  to  the  propellant  bum  rate  is 
during  the  actual  plasma  impact  and  there  are  no  secondary  bum  rate  augmentations.  This  was  shown 
by  Birk  [8]  who  also  did  some  experiments  with  chilled  propellants.  His  results  with  the  M30  propellant 
showed  that  the  bum  rate  at  -22  C  was  significantly  lower  than  those  samples  exposed  at  22  C.  Even 
though  he  didn’t  have  any  comparative  data  for  the  JA2  propellant,  one  can  conclude  that  a  similar 
effect  should  take  place  with  the  JA2  propellant. 

If  this  temperature  effect  had  a  considerable  impact  on  the  bum  rate  of  a  propellant,  then  this 
effect  must  be  quantified.  The  artillery  officer  in  charge  of  the  weapon  and  fire  guidance  must  know 
how  the  environment  is  going  to  impact  the  range  and  the  exterior  ballistics.  If  the  weapon  is  deployed 
into  an  environment  similar  to  Siberia  or  to  an  environment  like  the  Sahara  Desert  the  artillery  officer 
may  not  have  time  to  calibrate  the  devices  using  a  trial  and  error  method. 
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Information  about  the  firing  of  the  propellant  should  be  available  in  order  to  determine  if  the 
environment  will  affect  the  firing  and  to  what  extent  the  environment  will  impact  the  propellant. 

For  this  reason,  it  is  proposed  to  modify  the  bum  rate  equation  to  include  propellant  bed 
temperature  dependence  in  the  form 

BR(P,T)  =  aPbf(T ) 

All  of  the  previous  experiments  to  determine  the  constant  a  and  exponent  b  were  done  at  an  ambient 
temperature  assumed  to  be  20  °C,  but  the  propellant  temperature  was  not  monitored.  This  work  will 
determine  if  the  propellant  temperature  does  have  a  significant  impact  on  the  bum  rate  of  plasma  ignited 
propellant  and  if  so,  what  the  pressure-temperature  dependent  bum  rate  would  be.  It  is  assumed  that 
the  bum  rate  dependence  on  pressure  will  remain  as  a  power  law  but  it  is  possible  that  the  temperature 
dependence  could  follow  a  power  or  an  exponential  form 


BR  =  aPbTc 
or 

BR  =  aPb  exp (cT) 


The  bed  temperature  may  also  be  expressed  as  a  ratio  between  the  raised  temperature  and  the 
temperature  at  ambient: 


BR  =  aPb 


r  j\ 


\To  j 


where  T0is  the  ambient  temperature. 

The  expectation  is  that  there  will  be  propellant  temperature  sensitivity  in  the  bum  rate.  Beyer  [7] 
already  showed  some  effect  of  cooling  of  M30  propellant  and  this  side  note  found  in  his  work  must  be 
expanded  on.  However,  it  makes  sense  logically  that  the  temperature  of  the  propellant  bed  would  have 
an  effect  on  bum  rate  because  as  the  temperature  changes  in  the  bed,  so  does  the  energy  of  the  bed. 
Molecules  that  are  excited  are  more  likely  to  dissociate  from  each  other  when  ignited,  thus  increasing  the 
bum  rate. 

II.  Plasma  Flow-Field  Characterization 

Plasma  flow-field  characterization  and  distribution  of  plasma  parameters,  both  spatial  and  temporal, 
gives  an  accurate  picture  on  the  environment  where  the  propellant  bed  is  placed  within  during  exposure 
to  electrothermal  plasma.  With  controlled  bed’s  temperature  to  evaluate  the  sensitivity  of  the  propellant 
to  the  plasma  jet  at  elevated  temperatures,  it  is  necessary  to  determine  the  plasma  parameters  in  the 
plasma-flow  field.  Plasma  pressure  and  temperature  of  the  plasma  jet  outside  the  electrothermal  source 
must  be  known  by  direct  measurements.  These  measurements  can  help  determining  the  momentum  and 
energy  terms  at  any  chosen  location,  and  thus  the  plasma  impact  on  the  propellant  bed  can  be  expressed 
in  terms  of  the  plasma  parameters  at  chosen  location.  Before  conducting  temperature  sensitivity  studies, 
the  plasma-flow-field  experiments  were  conducted  and  the  pressure,  temperature,  density  and  velocity 
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distribution  in  the  flow-field  were  measured.  A  compact  multi-sensor  probe  has  been  designed  and 
constructed  to  measure  spatial  and  temporal  plasma  parameters.  These  parameters  are  the  total  and 
Static  plasma  pressure,  plasma  temperature,  plasma  number  density  and  plasma  Velocity.  All 
measurements  were  conducted  at  atmospheric  pressure  by  fusing  the  ET  source  with  aluminum  fuses. 
The  ET  plasma  is  thus  jetting  into  open  air  inside  of  a  cubical  chamber.  The  plasma  jet  expands  out  of 
the  ET  source  with  a  measured  expansion  angle  of  approximately  23°  from  the  ET  source’s  axis. 

II.  1  The  Compact  Multi-Sensor  Probe: 

The  compact  multi-sensor  probe  [23],  as  shown  in  the  Fig.  1  below,  has  two  absolute  pressure 
transducers  (Kistler  Model  60 1 B)  and  three  fiber  optic  cables.  The  sensor  head  is  arranged  such  that  the 
two  pressure  transducers  are  close  to  each  other  with  one  pointing  towards  incoming  plasma  flow  to 
measure  total  pressure,  and  one  situated  upright  to  measure  the  static  pressure. 


■T  Plasma 
Source 


/ 

-  V 

Compact  probe  head: 

2  pressure  transducers 
and  3  fiberoptics 
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Fig.  1  Schematic  of  the  concept  of  multi-sensor  probe  equipped  with  two  pressure  transducers  and 
three  fiber  optic  cables.  Positioning  of  the  pressure  transducers  is  shown  at  right. 


To  calibrate  the  pressure  transducers,  a  shock  tube  was  designed,  constructed  and  used  for  calibration. 
It  consists  of  a  high-pressure  chamber  and  a  barrel  section,  separated  by  a  burst  disc.  The  high-pressure 
chamber  is  connected  to  a  compressed  nitrogen  cylinder  and  set  to  the  desired  pressure.  The  valve  is 
shut  and  the  device  is  triggered  abruptly  via  a  plunger,  which  in  turn  forces  the  burst  disk  to  burst  and 
thus  forwarding  a  shock  wave  down  the  test  barrel.  Pressure  signals  are  recorder  on  a  LeCroy  9310 
digital  storage  oscilloscope.  It  was  determined  that  the  relation  between  the  true  pressure  versus 
measured  pressure  is  linear  and  is  expressed  by  PtrUc  =  1.03  ?measaKd  -  0.46.  The  three  fiber  optic  cables 
are  arranged  as  a  bundle,  with  one  fiber  optic  interfaced  to  an  optical  multichannel  analyzer  for  time- 
integrated  optical  emission  spectra  measurements,  and  two  other  fiber-optics  each  is  interfaced  to  a 
monochromator  and  a  photomultiplier  for  time-resolved  spectral  measurements.  Monochromators  were 
adjusted  to  521.82  nm  and  570.02  nm  wavelengths  (copper  lines).  Optical  emission  spectroscopy  [24- 
26]  was  used  to  measure  the  plasma  temperature,  assuming  plasma  is  at  LTE,  and  plasma-density  from 
optical  lines  broadening.  Three  axial  locations  were  chosen  for  measurements,  2.4,  3.4  and  4.4  inch 
from  source  exit.  At  each  axial  location,  three  radial  locations  are  defined  to  measure  distributions 
within  the  diverging  plasma  jet.  Figure  2  bellow  shows  the  source  exit  and  the  locations  of  compact 
multi-sensor  probe  measuring  positions. 
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Fig.  2  Schematic  showing  selected  locations  to  obtain  spatial  distribution  of  electrothermal  plasma 
parameters  down  stream  from  the  source  exit  and  within  a  20°  angle. 

Selection  of  the  measuring  locations  is  based  on  the  plasma  jet  divergence  as  the  plasma  leaves  the 
source.  The  angle  was  measured  by  placing  a  witness  plate  in  the  plasma  stream  along  the  axis  of  the 
device.  The  nine  selected  location  covers  the  solid  angle  of  divergence  and  provides  a  detailed 
information  on  the  spatial  distribution  of  the  plasma  jet  parameters. 

II.2  Plasma-Flow-Field  Distribution: 

Plasma  parameters  at  2.4,  3.4  and  4.4  inches  from  the  source  exit  were  measured,  and  at  three  radial 
locations  for  each  axial  position.  Figures  3  and  4  show  the  pressure  distribution  from  the  source  exit. 
Figure  3  shows  the  results  of  the  peak  total  pressure,  where  it  is  obvious  that  the  axial  profile  is 
approximately  flat  for  all  radial  locations,  while  the  radial  profiles  are  of  a  top-hat  distribution.  The 
measurements  imply  that  for  axial  locations  between  2  and  4  inches,  the  plasma  jet  undergoes  an 
essentially  isentropic  expansion,  while  the  edge  of  the  jet  is  subjected  to  mixing  process  with  the 
surrounding  gas. 

Figure  4  shows  the  distribution  of  the  peak  static  pressure,  where  the  general  trend  indicates  a 
decrease  in  the  static  pressure  as  the  jet  expands.  Of  course  the  largest  static  pressure  is  at  the  location 
closest  to  the  source  exit. 

Of  importance  is  the  distribution  of  the  plasma  temperature  and  density.  The  spatial  distribution 
of  the  plasma  temperature  is  shown  in  Fig.  5.  These  temperatures  were  obtained  from  time-integrated 
optical  emission  spectroscopy  using  copper  lines.  The  plasma  temperature  drops  quickly  from  2.1  eV 
right  at  the  source  exit  (axial  distance  =  0)  to  approximately  0.5  eV  at  2  inches  down  stream.  Further 
into  the  stream,  one  may  observe  that  the  plasma  jet  is  approximately  isothermal.  With  additional  data 
from  time-resolved  spectral  measurements  of  copper  lines,  it  was  obvious  that  the  plasma  temperature 
inside  the  jet  is  0.4  +/-  0.16  eV.  A  slight  increase  in  the  plasma  temperature  was  observed  at  4.4  inches 
from  the  source  exit.  Such  temperature  increase  is  attributed  to  the  expansion  of  the  plasma  jet. 


Fig.  3  Peak  total  pressure  in  the  expanding  plasma  jet 


Spatial  distribution  of  the  peak  static  pressure  (psig)  of  the  plasma  jet  as  measured  from  the  source  exit 


Fig.  4  Peak  static  pressure  in  the  expanding  plasma  jet 
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Fig.  5  Average  plasma  temperature  in  the  expanding  jet 


The  plasma  density  is  another  important  parameter.  It  was  calculated  from  the  peak  static  pressure 
and  the  time-integrated  plasma  temperature.  The  plasma  density  at  spatial  locations  of  measured 
pressure  and  temperature  is  shown  in  Fig.  6.  The  density  decreases  with  increased  axial  length  as  well  as 
increased  radial  position.  This  is  expected  for  an  expanding  plasma  jet  from  an  electrothermal  source.  It 
is  obvious  that  the  plasma  density  follows  the  static  pressure  distribution,  especially  that  the  temperature 
is  almost  same  inside  the  jet  and  thus  the  density  profile  follows  the  static  pressure  profile. 

The  average  speed  of  the  plasma  jet  was  obtained  from  the  plasma  arrival  time  to  the  sensors  at  the 
chosen  measuring  locations.  The  speed  of  the  jet  is  approximately  constant  and  close  to  1300  m/s  on 
axis.  The  slowing  of  the  plasma  near  the  jet  boundary  is  attributed  to  the  mixing  processes  with  the 
surrounding  air.  For  a  0.5  eV  ideal  gas,  the  sound  speed  is  3100  m/s,  and  thus  the  measured  speed 
indicates  that  the  jet  flow  is  sub-sonic  with  a  Mach  number  of  about  0.4.  The  time-resolved  spectral 
measurements  were  also  used  to  calculate  the  plasma  precursor,  from  which  the  precursor  plasma 
velocity  is  found  to  be  about  3000  m/s.  Fig.  7  shows  the  average  plasma  speed  in  the  expanding  jet. 

The  obtained  parameters  are  important  in  deciding  where  to  place  the  propellant  for  a  given 
experiment,  and  what  to  be  expected  in  terms  of  plasma  parameters,  at  any  selected  location. 


Fig.  6  Plasma  density  distribution  in  the  expanding  jet 


Fig.  7  Average  plasma  speed  in  the  expanding  jet 


III.  Temperature  Sensitivity  Studies 


10 


Temperature  sensitivity  of  the  propellant  bed  may  alter  the  bum  rates  under  plasma  injection.  When 
the  bed’s  temperature  is  elevated,  it  is  expected  that  the  bum  rates  would  increase  as  a  result  of 
increased  in-depth  heating  of  the  bed.  Also,  the  plasticizer  may  migrate  causing  alterations  in  the  bum 
rates.  A  temperature-control  unit  is  designed  such  that  the  bed’s  temperature  prior  to  plasma  injection  is 
maintained  constant.  The  bed  temperature  can  be  elevated  up  to  100°C  ,  however,  the  upper  limit  was 
set  to  50°C  [27], 

III.l  Experimental  Arrangement: 

A  circulating  cooling  unit  was  designed,  constructed  and  tested  for  operation.  Preliminary 
experiments  on  the  circulating  cooling  control  unit  have  shown  the  necessity  of  long  time  circulation  of 
the  fluid  in  order  to  maintain  temperature  stability.  A  modified  design  using  heating  elements  allowed  for 
a  faster  temperature  setting  prior  to  plasma  injection.  Figure  8  shows  the  experimental  arrangement 
using  heating  elements  to  fulfill  the  temperature  sensitivity  studies  with  elevated  temperatures  between 
20  to  50°C. 


Propellant 

sample 


Heating  element 


Thermocouple 


Heating 
element  leads 


Adjustable 
sample  holder 


Fig.  8  Schematic  of  the  experimental  setup  using  heating  a  heating  element  to  elevate  the  bed’s 
temperature  prior  to  plasma  injection  on  the  propellant  sample 


The  sample  holder  is  tube-shaped  and  has  a  seat  for  the  propellant  sample,  which  resides  on  the 
heating  element.  The  temperature  of  the  back  surface  of  the  propellant  is  monitored  via  a  thermocouple. 
The  sample  is  thin,  about  l-2mm  and  the  sample  temperature  is  kept  stable  to  the  desired  temperature 
via  the  heating  element  controller.  The  holder  has  a  fiberoptic  cable  to  view  the  plasma  at  the  plasma- 
propellant  interface.  The  entire  holder  is  adjustable  and  can  be  placed  at  any  location  on  the  axis  of  the 
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plasma  source.  Experiments  are  conducted  at  a  base  pressure  of  20  Torr  for  extinguished  bum 
(interrupted  bum)  such  that  the  bum  is  only  during  the  plasma  pulse  length  (100  pis  discharge).  The 
initial  and  final  weights  of  the  propellant  sample  are  recorded,  as  well  as  the  initial  and  final  weights  of 
the  plasma  source  liner  and  the  source  electrode.  For  each  shot,  the  discharge  current  is  fed  to  our 
TITAN  2-D  code  [28]  to  calculate  the  plasma  exit  parameters.  Plasma  temperature  at  the  plasma- 
propellant  interface  is  calculated  from  optical  emission  spectroscopy  measurements  using  copper  lines. 

A  block  diagram  of  the  experimental  arrangement  is  shown  in  Fig.  9,  where  PIPE  represents  the 
ETC  experimental  facility  powered  by  a  360pF,  lOkV  high  energy-density  capacitor.  The  block  diagram 
shows  experimental  components  and  data  acquisition  system,  thermocouples  and  fiber  optics  to  optical 
multichannel  analyzer. 


Fig.  9  Block  diagram  of  the  experimental  setup  for  temperature  sensitivity  measurements 


The  sample  holder  was  designed,  essentially,  to  be  a  rectangular  block  with  supporting  rods  that 
could  be  bolted  to  the  bottom  of  the  ETC  chamber  of  the  PIPE  facility.  The  heater  used  was  a  3/8" 
diameter  OMEGA  cylindrical  cartridge  heater  that  was  placed  on  a  shelf  and  held  in  place  by  an 
aluminum  plate.  A  small  copper  tube  was  attached  to  the  holder  by  epoxy  in  order  to  keep  the 
propellant  sample  from  sliding  against  the  pressure  and  provides  a  way  to  maintain  thermal  contact.  The 
initial  temperature  of  the  propellant  bed  was  monitored  by  a  thermocouple  pressed  against  the  exposed 
face  facing  the  plasma.  The  original  method  for  controlling  the  heater  output  was  an  Omega  Process 
Controller  that  controls  the  heater  temperature  using  an  output  relay.  This  controller  was  chosen 
because  it  would  keep  the  heater  temperature  essentially  constant,  and  could  be  monitored  and 
controlled  through  varying  the  voltage  across  the  heater.  However,  two  units  were  destroyed  in  a 
relatively  short  time  because  of  the  high-voltage  high-current  discharge.  A  large  potentiometer  was 
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added  to  the  system  and  used  to  power  the  heater.  The  potentiometer  regulates  the  voltage  on  the 
heater  at  a  given  time,  thus  controlling  the  current  and  temperature  of  the  heater.  This  setup  required 
careful  monitoring  of  the  rate  of  temperature  increase  so  that  the  experiment  would  be  performed  at  the 
time  corresponding  to  the  desired  propellant  bed  temperature.  It  was  possible  to  obtain  a  pseudo-steady 
state  temperature  with  this  device  after  an  hour  and  a  half  of  operation  because  the  heat  losses  to  the 
environment  would  eventually  balance  out  the  input  heat  flux  to  the  heater. 


III.2  Bed  Temperature  Sensitivity  Results 

Results  of  the  propellant  bed’s  temperature  sensitivity  have  shown  an  increasing  bum  rates  with 
increased  bed  temperature.  For  the  variation  of  bed  temperature  between  23  to  50°C,  the  best  fit  is 
showing  a  power  law  of  the  form:  Bum  rate  BR  =  A]  Tm,  where  the  bum  rate  is  measured  by  dividing 
the  interrupted  erosive  thickness  by  the  plasma  impact  time  (FWHM).  Here  in  this  equation,  T  is  the  bed 
temperature  in  °C,  A]  is  proportionality  constant  and  m  is  the  power  exponent.  This  equation  can  be 
introduced  to  the  bum  rate  equation  such  that  the  bum  rate,  including  bed  temperature,  may  be 
expressed  by:  BR  =  A  Pa  Tb,  where  the  first  term  (Pa)  represents  the  conventional  equation  for  the  bum 
rate  without  introducing  the  bed’s  temperature.  A  more  general  equation  may  be  introduced  as:  BR  =  A 
P  (  T/  Tambient) ,  and  thus  for  T  =  Talent ,  the  bum  rate  falls  to  the  known  equation.  Bearing  in  mind  that 
the  pressure  exponent  ‘a’  with  plasma  ignition  is  much  higher  than  its  value  for  conventional  ignition. 
The  following  graph  shows  the  experimental  results  for  elevated  temperature  regime  with  the  propellant 
situated  2  inches  from  the  source  exit.  When  the  propellant  is  closer  to  the  source  exit,  it  is  expected 
that  the  bum  rate  will  be  much  higher  because  of  the  pressure  effect  (momentum  effect)  on  the  bum 
rate. 


When  expressing  the  bum  rate  by  the  equation  Bum  Rate  BR  =  A  P“,  it  has  been  shown  for  JA-2 
that  the  bum  rate  in  conventional  ignition  is  given  by  BR  —  203.2x10 5  while  for  plasma  ignition  it 
was  shown  to  have  a  higher  exponent:  BR  =  8.865xl0'5  P1  376.  As  stated  above,  to  include  the  effect  of 
the  bed  temperature,  it  is  necessary  to  modify  the  bum  rate  equation  to  include  the  bed  temperature 
BR  (P,T)  =  a  Pb/T)  =  aPb  Tc,  where  T  is  the  bed  temperature.  The  bum  rate  may  be  expressed  in  terms 
of  normalized  bed  temperature  in  the  form 


BR  =  aPb 


\ToJ 


c 


Where  To  represents  the  ambient  (room)  temperature.  This  model  is  more  convenient  to  relate  the  bum 
rate  to  the  ratio  between  elevated  temperature  to  ambient  temperature,  or  could  be  expressed  as  the 
bum  rate  change  for  a  percentage  change  in  the  bed  temperature.  The  advantage  of  this  model  is  that  it 
directly  reduces  to  the  Saint  Roberts  bum  rate  at  ambient  temperature.  Figure  10  shows  experimental 
results  obtained  for  JA-2  bum  rates  versus  increased  bed  temperature  for  one  data  set  at  2-inches  from 
the  electrothermal  source  exit.  It  is  obvious  that  the  data  easily  fits  a  power  law  with  a  2.4767  exponent. 

The  bum  rate  data  for  all  selected  axial  locations  of  1,  2  and  2.75  inches  are  shown  in  Fig.  1 1  as  a 
function  of  the  normalized  temperature  (T/T0).  The  bum  rate  data  shows  that  the  1-inch  shots  have  the 
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highest  bum  rate  and  the  2-inch  the  lowest.  Initially  one  may  think  that  the  2.75-inch  data  should  be  the 
lowest  and  these  data  are  in  error.  However,  referring  to  the  pressure  data  obtained  by  the  compact 
multi-head  probe  as  previously  shown  in  section  II,  it  is  shown  that  total  plasma  pressure  actually 
increases  slightly  between  2  and  3  inches  as  a  result  of  momentum  transfer  and  expansion. 


Fig.  10  Bum  rate  of  JA-2  with  increased  bed  temperature.  Sample  placed  2  inches  from  the  exit  of  the 
ET  source.  Data  fit  to  power  law  with  2.4767  temperature  exponent 


Since  the  pressure  of  these  data  was  not  directly  measured,  a  method  of  evaluating  the  pressure  is 
to  compare  the  source  exit  pressure  with  previously  measured  pressures.  The  plasma  diagnostic  code 
TITAN  was  run  to  determine  the  barrel  exit  pressure  at  atmospheric  conditions  and  then  this  value  was 
fitted  to  the  values  obtained  by  compact  probe  measurements  (as  previously  described  in  section  II). 
TITAN  is  a  pseudo  2-D  code  that  calculates  the  plasma  parameters  at  the  source  exit  both  spatially  and 
as  a  function  of  time.  TITAN  was  also  used  with  base  pressure  of  20  Torr  (the  experimental  conditions) 
and  compared  to  compact  probe  results  to  see  if  the  chamber  pressure  had  any  significant  effect  on 
plasma  pressure.  TITAN  revealed  that  there  is  no  significant  change  between  the  plasma  expanding  at 
20  Torr  or  expanding  at  760  Torr  [6].  The  plasma  barrel  pressure  was  approximately  11,400  psi 
independent  of  the  initial  chamber  pressure.  This  value  was  fitted  to  compact  probe  data  at  2,3,  and  4 
inches  to  approximate  the  plasma  pressure  at  corresponding  distances.  The  pressure  values  are  shown  in 
Table  1  below. 


Table  1:  Plasma  Pressure  (Total) 


Distance  (inches) 

Pressure  (psia) 

1  1 

1403.6  1 

r  2 

295 

3 

304 
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Fig.  1 1  Bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions.  Normalized 
temperature  (T/T0)  is  used  in  this  graph 


Matlab’s  nonlinear  least  squares  function,  lsqnonlin,  was  used  to  determine  the  unknowns  in  the 
modified  bum  rate  equation  with  values  shown  in  Table  2. 
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Thus,  the  modified  bum  rate  equation  may  be  written  as: 


(  T  V  36 


BR  =  0.3 5 5P 


0.528 


\T.J 


JA2 


This  equation  fits  well  the  data  set  obtained  at  all  measured  locations,  as  seen  in  Fig.  12. 


Nomwdctd  TtAptnluri  TfT^ 


Fig.  12  Bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions.  Normalized 
temperature  (T/T0)  is  used  in  this  graph.  Data  fitting  used  normalized  bum  rate  equation. 

However,  it  can  be  seen  that  the  2-inch  data  is  not  well  fitted.  This  indicates  there  is  another  factor 
that  is  not  taken  into  account  in  the  above  equation.  The  equation  should  reduce  to  the  numbers 
obtained  by  Zaghloul  et  al  [5]  when  the  propellant  bed  reaches  ambient  temperature.  This  was  not  the 
case.  The  bum  rates  obtained  in  this  experiment  at  propellant  ambient  temperatures  were  much  higher 
than  those  obtained  previously  at  similar  pressures.  The  difference  is  attributed  to  the  way  the  nonlinear 
solver  fits  the  function.  A  quick  calculation  at  1 1 ,400  psi  and  propellant  bed  at  ambient  still  shows  an 
increased  bum  rate  with  plasma  injection  over  conventional  means.  With  plasma  injection  the  bum  rate 
at  ambient  temperature  is  49.2  cm/sec  verses  8.2  cm/sec  for  conventional  ignition.  The  equation 
obtained  by  Zaghloul  et  al  [5]  yields  a  bum  rate  of  33.72  cm/sec  with  plasma  injection.  In  general,  the 
new  modified  form  follows  the  same  trend  even  though  the  pressure  exponent  is  not  as  large. 


This  is  due  to  the  large  difference  in  the  magnitude  of  the  constant  ‘a’  and  the  pressure  range  looked 
at.  The  equation  in  this  form  does  provide  a  good  application  oriented  equation  because  impact  pressure 
is  a  common  measurable  value  in  cannon  ballistics.  It  is  also  possible  to  remove  the  pressure  dependence 
by  taking  the  ratios  of  the  bum  rates  at  a  raised  bed  temperature  to  the  bum  rate  at  ambient  bed 
temperature.  Thus,  the  equation  may  then  be  written  as: 

( br \  _  aF^r  _  J ’  r_Y 

(BR\  “  aP'TJ  =  \jo)ja 

where  C  is  a  constant.  This  equation  is  valid  for  other  ignition  means,  e.g.  conventional  or  laser, 
because  it  depends  only  on  the  propellant  bed  temperature. 

Each  bum  rate  data  set  was  analyzed  and  was  normalized  to  the  ambient  propellant  bum  rate.  The 
values  used  for  the  normalization  were  the  bum  rates  at  T=20°C  or  the  average  of  the  coolest  ambient 
data  point  available.  The  1-inch  data  was  normalized  to  34.325,  2  inches  to  8.566,  and  2.75  inches  to 
3.485  cm/sec.  Normalized  bum  rate  is  shown  in  Fig.  13. 


Fig.  1 3  Normalized  bum  rate  of  JA-2  with  increased  bed  temperature  for  all  samples’  positions. 
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It  can  be  seen  from  Fig.  13  that  the  normalized  bum  rates  at  1  and  2  inches  are  almost  the  same 
while  the  normalized  bum  rate  at  2.75  inches  increases  significantly. 


The  normalized  equation  was  transformed  into  a  log  form  and  a  linear  fit  was  performed. 


In 


f  (M),  ' 


(BR) 


T,  J 


=  ln(C)  +  yln 


f  y  ^ 

V  ^0  Jja2 


The  fit  parameters  are  shown  in  Table  3  for  1, 2,  and  2.75  inches  from  the  electrothermal  source  exit. 


Table  3:  Temperature  Exponent  y  and  constant  C  for  1,  2  and  2.75  inches 


C 

Y 

1-inch 

0.6110 

1.9721 

2-inches 

0.5638  ~ 

2.4767 

2.75-inches 

1.0352 

3.3312 

These  fits  are  shown  in  Figures  14,  15  and  16. 


Fig.  14  Normalized  Bum  Rate  Fit  for  d=l  inch  from  the  source  exit 
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Fig.  15  Normalized  Bum  Rate  Fit  for  d=2  inch  from  the  source  exit 


Fig.  16  Normalized  Bum  Rate  Fit  for  d=2.75  inch  from  the  source  exit 


These  figures  show  that  the  propellant  bed  temperature  augments  the  data  more  the  farther  it 
extends  beyond  the  barrel,  which  indicates  the  bum  rate  is  also  dependent  on  the  distance  from  the 
source.  The  conclusion  here  is  that  the  plasma  parameters  do  play  an  important  role  in  propellant  bum 
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and  need  to  be  incorporated  into  the  equations.  An  additional  functional  form,  with  dependence  on  bed 
temperature  and  distance  from  the  source  section  is  proposed. 


BR 

BR0 


f 


=  C 


\To  JjAl 


(  J  V 


\doJ 


where  the  naught  values  are  taken  at  an  arbitrary  distance  from  the  source. 

Since  P  =  I  nkTp,  the  pressure  may  then  be  substituted  in  the  bum  rate  equation  to  include 
plasma  parameters  in  lieu  of  the  pressure: 


BR  =  a(nkT '  ) 


»{l}r 


where  n  is  the  plasma  density,  k  is  the  Boltzmann  constant,  and  Tp  is  the  plasma  temperature. 

Plasma  temperature  can  be  related  to  plasma  heat  flux  utilizing  the  Blackbody  assumption  such  that 
the  plasma  temperature  could  be  expressed  in  terms  of  the  blackbody  heat  flux  q" : 


T  = 


Vy 

UcJ 


where  e  is  the  emissivity  and  a  is  the  Stefan-Boltzmann  constant.  For  simplicity,  the  emissivity  s  is 
assumed  to  be  equal  to  1 ,  and  thus,  substituting  for  the  plasma  temperature  yields  an  expression  for  the 
bum  rate  in  terms  of  radiation  heat  flux: 


where  A  is  a  constant. 


BR  =  Anb 


Optical  emission  spectroscopy  has  been  performed  to  determine  the  plasma  temperature  and  density. 
Copper  lines  were  used  to  obtain  plasma  temperature  and  number  density.  The  lines  used  for  the 
temperature  analysis  and  their  parameters  are  shown  in  Table  4. 


Table  4:  Copper  lines  used  for  Spectroscopic  Analysis 


Spectrum 

Wavelength  (nm) 

Ei(eV) 

gA(10V) 

Cu  I 

510.554 

3.82 

0.080 

Cu  I 

515.324 

6.19 

2.4 

Cu  I 

521.820 

6.19 

4.5 

Cu  I 

529.250 

7.74 

0.436 

Cu  I 

570.020 

3.82 

0.0096 

Cu  I 

578.213 

3.79 

0.033 

InttnfiV  (Hu.)  lnitnfity(M.)  tahroilyCM.) 
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Three  spectral  samples  are  shown  for  each  distance  in  Figures  17,  18  and  19,  to  demonstrate 

reproducibility  of  the  spectral  data,  thus  indicating  similar  properties. 


X  104 


(n«) 


Fig.  17  Sample  three  d  =  1-inch  shot  Spectra 
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The  PEAKFIT  [29]  program  was  used  to  determine  all  measured  peaks  in  the  spectra.  It  can 
determine  all  the  peaks,  obvious  and  hidden,  and  outputs  them  and  their  areas  to  file  for  analysis.  Table 
5  shows  plasma  temperature,  density  and  the  corresponding  radiation  heat  flux  at  the  three  locations  1, 
2  and  2.75  inches  from  the  source  exit. 


Table5:  Average  Plasma  Temperature,  Heat  Flux,  and  Density 


Distance 

_ (in) 

Temperature 

_ (eV) 

Heat  Flux  *108 
W/m2 

Density  *1024 

m'3 

1  1 

0.846 

5.237 

3.575 

2 

0.687 

2.289 

3.406 

2.75 

0.821 

4.666 

2.614 

The  purpose  for  determining  these  temperatures  is  to  attempt  to  determine  the  bum  rate 
dependence  on  plasma  parameters  as  well  as  propellant  bed  temperature.  A  bum  rate  ratio  can  be 
determined  between  the  distances  that  will  allow  correlation  with  number  density  and  plasma  heat  flux. 
The  ratio  can  be  expressed  as: 


(BR) 
( BR\ 
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'  n  ' 
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JA2 


The  only  question  about  this  equation  is  what  the  normalizing  condition  should  be.  The  reasonable 
choice  is  the  plasma  density  because  it  dominates  the  plasma  temperature  in  the  kinetic  pressure  term  P 
=  I  nkT.  For  this  reason,  all  naught  values  are  those  of  the  2.75  inch  data  at  ambient  temperature  so 
that  the  bum  rate  ratio  would  reduce  to  one  at  that  point.  Utilizing  the  values  obtained  for  plasma 
temperature  and  density  allows  the  unknowns  of  the  above  equation  to  be  determined  using  Matlab’s 
nonlinsq  function  and  are  shown  in  Table  6. 

Table  6:  Parameters  of  the  modified  equation  with  dependence  on  plasma  radiation  heat  flux  and  number 


density,  and  the  normalized  bed  temperature  with  respect  to  ambient 


Parameter 

A . . *  _ _  _ 

Value 

A 

1.8315 

b 

2.6114 

_ 1 _ 

2.8365 

Hence,  the  bum  rate  equation  can  now  be  written  as: 


m 

(BR)a 


1.8313 


\  2.6114/  n  0.6529  / 

a  t 


Vno  J 


Wo’ 


yT0. 


2.8365 

JA2 


This  Equation  provides  an  insight  into  the  effect  of  plasma  density  and  radiation  heat  flux  on  the 
bum  rate  of  the  propellant  bed.  It  is  now  possible  to  predict  bum  rate  augmentations  based  on  the 
plasma  density  and  heat  flux  variations  as  well  as  propellant  bed  temperature. 


IV.  Arc  Channels 
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IV.  1  Arc  Channel  Experiment 

Experimental  arrangement  to  initiate  an  arc  channel  in  a  propellant  sample  was  arranged  using  a  high 
voltage  pulsing  device.  Such  arcs  generated  by  HV  pulsers  are  usually  easy  to  generate  at  atmospheric 
conditions.  The  arc  may  take  any  shape  from  one  pulse  to  another  depending  on  the  triggering  circuitry 
and  the  atmospheric  conditions  during  the  time  of  measurements.  However,  with  short  arcs  generated 
between  thin  electrodes  at  small  separating  distance,  consistency  of  the  arc  parameters  could  be 
obtained.  Figure  20  shows  the  experimental  arrangement  of  developing  an  arc  channel  in  a  propellant 
sample.  The  sample  is  situated  between  two  electrodes  to  the  output  of  a  high  voltage  pulser.  The  pulser 
may  be  triggered  once  or  multiple  pulses  may  be  produced  at  a  selected  repetition  rate. 


insulated  leads 


Trigger 

Pulse 


/ 

HV  pulser 

V 

Propellant 

sample 


Pulser 

Conductive  leads 


Fig.  20  Experimental  arrangement  to  develop  an  arc  channel  in  a  propellant  sample  at  atmospheric 
conditions.  Single  or  multiple  pulses  may  be  produced. 


The  propellant  sample  resides  between  the  pulser  leads  such  that  the  arc  is  initiated  through  the 
central  perforation  of  the  JA-2  sample.  This  will  initiate  the  arc  with  its  core  inside  the  propellant  as 
illustrated  in  Fig.  21(A),  where  the  arc  core  joule-heats  the  inner  wall  of  the  propellant.  Also,  an 
exploding  wire  may  be  inserted  in  the  central  perforation,  Fig.  21(B)  and  connected  to  the  pulser  to 
initiate  the  arc  via  wire  explosion. 


energetic  material 
‘Propellant  bed" 


arc  initiation  exploding  wire 


A  B 

Fig.  21  Arc  initiation  inside  the  central  perforation  (A),  and  by  fuse-operated  arc  (B)  inside  the  central 
perforation  of  a  propellant  sample. 


The  arc  voltage  can  be  increased  to  up  to  30kV  with  arc  current  up  to  70 A,  and  pulse  length  of 
approximately  10-20  ns,  giving  rise  to  arc  energies  of  several  of  mJ.  Typical  arc  current  has  a  peak  of 
about  1.0  A  at  a  peak  discharge  voltage  of  12kV  for  most  experiments.  The  integrated  energy  of  a 
typical  arc  is  less  than  120  mJ,  which  makes  it  perfect  to  simulate  the  mJ  energy  range  of  arc  channels 
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developing  in  a  propellant  bed.  Controlling  the  arc  energy  in  the  present  pulser  is  via  the  arc  voltage 
which  provides  a  wide  range  of  arc  energy  between  few  milli-joules  to  several  joules.  The  arc  current  is 
sharp  enough  to  closely  simulate  initiation  of  narrow  arc  channels  at  considerably  short  period  of  time. 
The  arc  current  is  measured  by  a  Pearson  coil  and  arc  voltage  is  measured  by  a  compensated 
capacitively-coupled  high  voltage  probe.  Figure  22  shows  a  schematic  of  the  experimental  arrangement 
tor  milh-joule  to  several  joules  arc  discharge. 


Fig.  22  Schematic  of  experimental  arrangement  for  low-energy  arc  channels  in  solid  propellant  sample. 

Experimental  data  have  shown  a  weight  change  of  ±0.03  mg  for  such  low-energy  arcs  which 
prompts  the  need  to  either  increase  the  arc  energy  or  to  pre-stress  the  sample  by  compression  The 
choice  of  mcreasmg  arc  energy  has  been  chosen  to  investigate  effect  of  arc  energy  on  initiation 
threshold.  Arc  energy  between  260  and  1050  joules  were  used  and  the  corresponding  mass  loss  was 
measured.  Figure  23  shows  experimental  mass  loss  data  versus  arc  energy,  and  a  linear  fit  to 
experimental  data.  The  linear  fit  Am  =  2  x  10-4  (E  +  282),  where  Am  is  the  mass  loss  in  mg  and  E  is  the 
arc  energy  m  joules.  The  linear  fit  has  R2  =  0.8917,  almost  0.9.  This  medium-range  arc  energy  shows  an 
increased  erosive  bum  with  increased  energy  of  the  arc,  and  it  is  consistent  with  previous  results  of 
plasma  enhanced  bum  rates  of  same  propellant  [4,5]. 


Mw»  Lou  vo.  Input  Enorgy 


Fig.  23  Measured  weight  loss  of  propellant’s  sample  versus  arc  energy  showing  increased  bum  rate  with 
mcreased  energy.  The  linear  fit  has  R2  ~  0.9. 


IV.2  Arc  Channel  Model 
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Arc  channel  may  be  modeled  as  a  confined  capillary  discharge,  which  is  similar  to  modeling 
electrothermal  plasmas  in  ablation-controlled  arc  regime  [1 1,14,28]. 

The  electron-ion  collision  frequency  is  necessary  to  compute  e-i  momentum  transfer.  The 
Coulomb  logarithm  in  such  arc  regime,  with  high  density  and  low  temperature,  needs  be  modified  and  the 

approximation  model  proposed  by  Zollweg-Liebermann  has  formula  [30]  In  A  =  In  ^1  +  1 .4A2m  ,  where 
the  term  Am  is  a  modified  Coulomb  parameter,  and  is  equal  to: 

_  ^XD+X) 


Here  AD  is  the  Debye  length  (m),  A+  is  the  mean  ionic  radius  (m),  and  bQ  is  the  average  thermal  impact 
parameter  (m)  for  90°  scattering  of  incident  electrons.  The  mean  ionic  radius  is  given  by: 


'max 


This  summation  is  taken  over  all  possible  ionization  states  (/?,  is  the  number  density  of  ionization  state  i ), 
including  neutrals  (represented  by  z  =  1).  Consequently,  the  electron-ion  collision  frequency,  including  the 
modified  Coulomb  logarithm,  will  be  given  by  [28,30]: 

_  _  38 Zeffnee2  In  ^1  +  1.4A2m 
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with  Zeff  is  the  effective  charge  state  of  the  plasma,  ye  is  the  correction  factor  for  e-e  collisions,  and  me  is 
the  rest  mass  of  the  electron.  The  conductivity  of  the  plasma  ‘arc’  is  finally  given  by: 
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Here  Qe0  is  the  momentum  transfer  cross  section  for  electron-neutral  collisions.  The  conductivity  directly 
influences  the  volumetric  heating  of  the  plasma,  which  is  assumed  to  be  composed  only  of  Joule  heating: 

j2 

q"'=  — 

J  is  the  magnitude  of  the  vector  quantity  of  the  current  density  (A/m2),  and  a  is  the  plasma  conductivity 
(S/m).  This  volumetric  heat  generation  rate  (W/m3)  then  goes  into  the  energy  equation,  which  is  coupled 
to  the  momentum  and  continuity  equations.  They  are  given  below,  in  the  order  of  energy,  momentum, 
and  then  continuity: 


dt 

d(pv) 


+  V-(puv)=-PV-v+V-AVT  +  q’" 
+  V-(pvv)  =  -VP  +  V-<rz+Jx5 


dt 
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~^’(pv)  =  o 

A  number  of  variables  are  introduced  here.  Thermodynamic  variables  include  the  pressure  (Pa), 
temperature  (K),  specific  internal  energy  (J/kg),  and  physical  density  (kg/m3),  which  are  P,  T,  u,  and  p, 
respectively.  Flow  variables  include  the  velocity  (m/s),  conductivity  (sum  of  thermal  and  radiative 
conductivities,  W/m  •  K),  and  shear  stress  (N/m2),  given  as  v ,  A,  and  c?z ,  respectively.  Finally,  the 

electromagnetic  terms  are  the  current  density  (A/m2)  and  magnetic  field  strength  (T),  given  as  J  and 

B ,  respectively  [28].  The  code  results,  for  an  assumed  flat-top  arc  current  over  a  50  ps  period,  inside  a 
polycarbonate  tube  that  simulates  central  perforation  of  a  propellant  have  shown  higher  temperatures 
than  observed  ones.  Figure  24  shows  the  temporal  plasma  arc  temperature  as  calculated  by  the  code.  It 
is  suggested  that  the  approximation  of  the  Coulomb  logarithm  is  not  accurate  and  results  in  over¬ 
estimation  of  the  plasma  arc  temperature.  The  code  results  for  aluminum  capillary,  simulating  fuse 
operation,  is  shown  in  Fig.  25,  where  the  temperature  peaks  to  almost  same  value  as  that  of 
polycarbonate.  The  difference  is  in  the  fester  temperature  decay  for  the  aluminum  case,  due  to  heat 
conduction  to  the  channel  wall.  However,  it  is  over-predicting  the  temperature.  A  modification  to  the 
Coulomb  logarithm  should  be  introduced  to  account  for  plasma  non-ideality  with  better  approximation 
or  a  complete  non-approximated  solution.  This  modification  will  be  investigated  in  a  future  phase. 
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Fig.  24  Calculated  plasma  arc  temperature,  in  a  polycarbonate  capillary  tube,  as  a  function  of  time  for  a 
flat-top  current  profile. 
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Fig.  25  Calculated  plasma  arc  temperature  for  an  aluminum  capillary,  simulating  aluminum  fuse,  as  a 
function  of  time  for  a  flat-top  current  profile. 


V.  Conclusion 


The  plasma  flow-field  parameters  were  measured  showing  the  distribution  of  essential  parameters, 
such  as  pressure,  temperature,  number  density  and  average  speed.  The  importance  of  such  plasma- flow- 
field  characterization  is  to  help  in  the  design  of  electrothermal-chemical  devices  by  knowing  the 
distribution  of  the  plasma  jet  parameters.  c 


It  is  concluded  that  there  is  a  power  law  bum  rate  dependence  on  propellant  bed’s  temperature.  This 
has  been  shown  experimentally  and  has  a  reasonable  physical  meaning  because  of  the  extra  energy  given 
to  the  propellant  bed  due  to  the  heat  energy  input.  The  functional  form  of  bum  rate  equation  haa  been 
expressed  in  various  forms  but  always  includes  the  bed’s  temperature.  A  useful  form  is  that  of  the 
normalized  bed  temperature  with  respect  to  ambient  temperature  and  a  dependence  on  the  pressure. 
Another  form  is  the  ratio  form  between  bum  rates  at  elevated  temperature  to  that  at  ambient.  The  form 
that  includes  plasma  parameters  in  terms  either  of  plasma  temperature  or  radiation  heat  flux  is  of  more 
importance  in  order  to  characterize  effect  of  plasma  augmentation  mechanism  on  bum  rates. 

Arc  modeling,  using  approximation  to  the  Coulomb  logarithm,  has  shjown  over-prediction  of  the 
plasma  arc  temperature.  This  over-prediction  has  been  proven  for  the  case  of  polycarbonate  and 
aluminum.  The  approximation  affects  the  calculation  of  the  electron-ion  collision  frequency  and 
consequently  the  plasma  arc  conductivity,  which  directly  affects  the  joule  heating  term  in  the  energy 
equation.  A  better  approximation  of  the  modified  Coulomb  logarithm,  or  a  complete  non-approximated 
form  should  be  introduced  to  account  for  plasma  non-ideality,  which  will  be  investigated  in  a  future 
phase. 
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